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We examine the Big Bang nucleosynthesis (BBN) process in the presence of small-scale baryon 
inhomogeneities. Primordial abundance yields for D, 4 He, 6 Li, 7 Li, 9 Be, and n B are computed 
for wide ranges of parameters characterizing the inhomogeneities taking account of all relevant 
diffusive and hydrodynamic processes. These calculations may be of interest due to (a) recent 
observations of the anisotropies in the cosmic microwave background radiation favoring slightly 
larger baryonic contribution to the critical density, Ob, than allowed by a standard BBN scenario and 
(b) new observational determinations of 6 Li and 9 Be in metal-poor halo stars. We find considerable 
parameter space in which production of D and 4 He is in agreement with observational constraints 
even for Qf,h 2 a factor 2-3 larger than the fif, inferred from standard BBN. Nevertheless, in this 
parameter space synthesis of 7 Li in excess of the inferred 7 Li abundance on the Spite plateau results. 
Production of 6 Li, 9 Be, and 11 B in inhomogeneous BBN scenarios is still typically well below the 
abundance of these isotopes observed in the most metal-poor stars to date thus neither confirming 
nor rejecting inhomogeneous BBN. In an appendix we summarize results of a reevaluation of baryon 
diffusion constants entering inhomogeneous BBN calculations. 



I. INTRODUCTION 

The possibility that cosmic baryon number fluctua- 
tions may have existed on small scales in the early uni- 
verse has received considerable attention between the 
late eighties and mid nineties Such fluctuations in 

baryon number would have impact on the production of 
light elements during BBN provided the baryonic mass 
of individual lumps exceeds Mb > 10~ 21 Mq . It was 
speculated that production of inhomogeneities could re- 
sult during a first-order QCD phase transition or even 
possibly during a scenario of electroweak baryogenesis. 
Initially it was hoped for that such scenarios could make 
BBN consistent with — 1 and therefore eliminate the 
need for "exotic" non-baryonic dark matter. Detailed cal- 
culations revealed that inhomogeneous BBN (hereafter, 
IBBN) scenarios may not be consistent with a universe 
closed in baryons due to considerable overproduction of 
7 Li and/or 4 He. At present, = 1 seems also hardly 
desirable because of a variety of other arguments, such as 
the cosmological baryon budget |(| , and the success of a 
structure formation scenario employing cold dark matter, 
among others. 

Recently, observations of the cosmic microwave back- 
ground radiation (hereafter; CMBR) on intermediate an- 
gular scales by the BOOMERANG and MAXIMA bal- 
loon missions have achieved unprecedented accuracy . 
These missions have allowed for a first stab at an es- 
timate of a number of cosmological parameters such as 
the total cosmic density parameter fi tot and Qb, among 
others. Common to both studies is the observation of a 
relatively suppressed CMBR power spectrum on scales 
where a secondary peak in the spectrum is anticipated 
compared to the power at the location of the first peak. 
Though preliminary, the conclusion of a number of au- 



thors is that, within the parameters commonly allowed to 
be varied, an increased fib could most easily account for 
such a suppression ||. (for alternative explanations cf. 
to Ref. j9|). This has lead to the preliminary claim that 
fib as inferred from CMBR anisotropy observations may 
be in conflict with the best estimate flbh 2 0.02 ± 0.002 
from SBBN in particular, the CMBR data would 

prefer flbh 2 w 0.03 M (h is the Hubble constant in units 
of 100 km s _1 Mpc ). Note that even though, at first 
glance the deviation between these two values seems rela- 
tively small, it is clear that a baryonic density parameter 
of SI5 w 0.03 hr 2 can not be achieved within a SBBN 
scenario. For such large fib, SBBN production of deu- 
terium can neither account for the deuterium as observed 
in quasar absorption systems |R|[ll|] , nor for the inferred 
D abundance in the presolar nebula and only barely for 
the D as observed in the local interstellar medium. 

Newly developed high-resolution spectrographs (such 
as UVES on the VLT) allow for a significant increase 
in the number of stars with claimed detections for the 
elements 6 Li and 9 Be. Whereas for a long time there 
had been only two claimed 6 Li/ 7 Li Q| detections in low- 
metallicity PopII halo stars, this number is/will rapidly 
increase in the immediate future. Moreover, 6 Li/ 7 Li de- 
tections have now also been claimed for disk stars at 
relatively high metallicities The preliminary pic- 

ture which emerges is that ^i/H abundances in stars 
are remarkably similar over a wide range in metallici- 
ties, though interpretation of the data has to account for 
the possibility of stellar 6 Li astration. Recently, there 
has been an interesting 9 Be/H detection within the at- 
mosphere of a very low-metallicity star fl6|| . The 9 Be/H 
abundance in this star is higher than expected from ex- 
trapolation of the approximately linear 9 Be/H versus 
[Fe/H] relation such that this observation may represent 
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tentative evidence for a flattening of the 9 Be/H versus 
[Fe/H] slope at metallicities below [Fe/H] < —3. 

In light of the above, it seems worth reinvestigat- 
ing BBN with an inhomogeneous baryon distribution. 
Whereas production of 6 Li and 9 Be in standard BBN is 
essentially negligible, it is known that production of these 
isotopes in IBBN may be significantly enhanced Jl7| . 
Furthermore, it should be of interest not only to find 
out of how much the upper limit on flf, in IBBN may be 
relaxed compared to that from standard BBN, but also 
in how much of the parameter space, characterizing the 
inhomogeneities, IBBN abundance yields may agree with 
observational constraints. 



II. INHOMOGENEOUS BIG BANG 
NUCLEOSYNTHESIS CALCULATIONS 

We have performed detailed numerical computations of 
IBBN by employing the IBBN code described in Ref || . 
This code treats all the relevant baryon diffusion of neu- 
trons, protons, and lighter nuclei. In the Appendix we 
summarize the employed baryon diffusion constants for 
protons and neutrons, which includes a reevaluation of 
some diffusion constants and a correction for mistakes in 
the literature. The employed code is still the only exist- 
ing code with a detailed treatment of the effects of pho- 
ton diffusion and hydrodynamic expansion on the evo- 
lution of high-density regions j^n)). It is known that 
these dissipative processes operating at lower tempera- 
tures T < 30keV may affect the predicted abundances 
of 6 Li, 7 Li, 9 Be, and 11 B in some part of the param- 
eter space (in particular, for compact high-density re- 
gions). For example, 7 Li produced in form of 7 Be may 
be prematurely destroyed by the reaction sequence 7 Be 
(n,p) 7 Li (p,a)a when enough neutrons may be deliv- 
ered to the high-density regions where most of the 7 Be is 
produced. The magnitude of this process depends on the 
efficiency of hydrodynamic expansion which increases the 
surface area of high-density regions but also on the cor- 
rect neutron- and proton- diffusion constants at low tem- 
peratures. Note that the distribution (and diffusion) of 
protons affects the diffusion of neutrons through neutron- 
proton nuclear scattering Q . 

We have updated the nuclear reaction rates employed 
in the IBBN code from those based on the compilation 
by Caughlan & Fowler, as described in Smith et al. [5o| , 
to include the improved charged nuclei induced reactions 
as compiled by the NACRE collaboration [|TJ . Note that 
the modifications in predicted abundances, when the cen- 
tral values of the improved nuclear reaction rates of the 
NACRE compilation are employed, are fairly small for 
2 H and 4 He but can be in the ~ 20% - 30% range for 7 Li, 
6 Li, 9 Be, and n B (cf. to Ref. |§]). Additional uncer- 
tainties in the computed abundances arise from appre- 
ciable error bars quoted in the NACRE compilations for 
a few reactions, such as D(p,7) 3 He, 3 He(a,7) 7 Be, and 



D(a,7) 6 Li. In the context of SBBN, these additional 
uncertainties are of similar magnitude to those quoted 
above, with the exception of the 6 Li abundance which 
is subject to very large uncertainties of ~ factor 3-4 in 
either direction [^,^2| . Though it is beyond the scope of 
the present work to present a detailed systematic anal- 
ysis of uncertainties in the prediction of abundances in 
IBBN scenarios due to reaction rate uncertainties, we 
will comment below if such uncertainties could impact 
our main conclusions. Note that all 4 He mass fractions 
Y p presented below are corrected by AY p = +0.0049 to 
account for a variety of physical effects as detailed in 
Ref. |§. 

From the multitude of conceivable initial conditions for 
the baryon inhomogeneities (including stochastic ones as 
treated in Ref. j2j| ) we chose a regular lattice of spherical 
symmetric domains, approximating the possible outcome 
of baryon fluctuations generated during a first-order (e.g. 
QCD) phase transition around the shrinking bubbles of 
high-temperature phase. The spherical computation do- 
main is then characterized by its physical length, hoo, 
specified at temperature of T = 100 Me V (specifically, 
1 m at T = 100 MeV is to be understood as a length 



of 5.96 x 10 m at the present epoch 26 ). Within this 



domain we assume a region of high baryon density with 
baryon-to-photon ratio rjh occupying volume fraction fy 
and a low density region at r]\ = rjh/R occupying the 
remainder of the volume, with an initial discontinuity 
at the boundary of both regions (which softens after 
some baryon diffusion). This yields an average baryon- 
to-photon ratio 



v = fvRm + (i - fy) -qi 



(i) 



Given these initial conditions there exist still four pa- 
rameters to be specified, namely ij, ixooi fv, and R. The 
initial parameter space is reduced by assuming fyR = 
200. Physically fyR ^S> 1 corresponds to essentially all 
baryons residing in the high-density region and none in 
the low-density region, such that for fyR > 10 one ob- 
tains results essentially independent of the exact value 
of this parameter combination. Though the opposite 
limit, fyR -C 1, may be interesting for the production 
of significant amounts of isotopes with nucleon number 
A > 12 [^7j, in much of the parameter space it yields 
only minor changes in the D, 4 He, and 7 Li as compared 
to a SBBN scenario at the same rj. Our calculations 
employ two different initial "geometries": (a) spherical 
condensed - where the high-density region resides at the 
center of the spherical domain and (b) spherical shell 
- where the high-density region occupies a shell at the 
outer edge of the computational domain. These spherical 
domains are finite-differenced into 24 zones. By increas- 
ing the number of zones we estimate that the relative 
error in predicted abundances does not exceed ~ 0.5% 
for 4 He, ~ 3% for D, and ~ 10% for the shown isotopes 
with A > 7. 
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III. RESULTS AND DISCUSSION 

Figures 1-4 show computed abundance yields in IBBN 
scenarios for the isotopes of D, 4 He, 7 Li, as well as 
6 Li, 9 Be, and n B as a function of the length scale of 
the domains (approximately corresponding to the mean 
separation between fluctuations) for differing fib and a 
wide range of parameters describing the baryon inhomo- 
geneities. The choice of the parameter space for which 
abundance yields are shown is supposed to bracket most 
potentially interesting (taking values of Vl^h 2 = 0.012, 
0.025, 0.038 and 0.051 shown by the solid, dotted, short- 
dashed, and long-dashed lines in each figure, respec- 
tively), as well as to illustrate the general trends of 
changing "geometry" (Figures 1 and 2 are for spherical 
condensed fluctuations, whereas Figure 3 and 4 are for 
spherical shells) and changing volume fractions fy of the 
high-density regions. We have deliberately not indicated 
observationally inferred limits on the primordial abun- 
dances in these figures, as these are likely to change over 
the course of time, and since we are more interested in 
a qualitative understanding of IBBN abundance yields 
and their potential agreement/disagreement with obser- 
vationally inferred abundance limits. 

It is well known that there exists an "optimum" length 
scale where the 4 He yield may be lower than in a SBBN 
scenario at the same Qb due to the fact that neutrons 
may diffuse out of the high-density regions, subsequently 
decaying in the low-density region before they may be 
incorporated into 4 He. This effect is more pronounced 
when fy is small (Figures 1 and 3) since back-diffusion 
of neutrons into the high-density regions is less efficient. 
Similarly, one finds an "optimum" ^ioo where 7 Li pro- 
duction is minimized, generally somewhat smaller than 
the length scale for minimum 4 He production. Never- 
theless, even at this "optimum" distance the 7 Li yields 
are often higher, at best somewhat lower, than the 7 Li 
yields in a SBBN scenario at the same fib and they in- 
crease with decreasing fy pg}] . These trends are due to 
the two different production mechanisms for 7 Li (direct 
production in the low-density region and production of 
7 Be in the high-density region) and the relative efficiency 
of these mechanisms at either lower or higher r\ than the 
approximate r\ inferred from SBBN. 

There are currently two mutually inconsistent obser- 
vationally inferred values for the 4 He mass fraction, i.e. 
a high value Y p sw 0.244 |2{| and a low value Y p sa 
0.234 I^HH . It also becomes more and more appreciated 
that the inference of Y p from observations of HII regions 
is subject to systematic errors of possibly considerable 
magnitude. Much progress has been made in the determi- 
nation of primordial D abundances in quasar absorption 
line systems (Q ASs) . There are now several QASs seem- 



ingly indicating low D/H ps 2.5 — 4x 10~ 5 
one which favors high D/H « 2 x 10~ 4 
this, one should probably demand from a sucessful BBN 
scenario to have Y p < 0.25 (conservative) and low D/H 



|13,|11| and only 
|§. In light of 



« 2 — 5 x 10 5 . (Note that even within the context of 
a SBBN scenario, a D/H «3x 10 -5 abundance implies 
seemingly uncomfortably high Y p m 0.247 0,0. The 
figures illustrate that one may find considerable IBBN 
parameter space where these requirements on the pri- 
mordial 4 He and D/H abundances may be met, even for 
flbh 2 as large as ~ 0.05. This is due to IBBN scenarios 
often yielding less 4 He and more D production than a 
SBBN scenario at the same fib- 

Nevertheless, these considerations disregard observa- 
tional limits on the 7 Li/H abundance. Typical 7 Li yields 
in the IBBN parameter space which agree with obser- 
vational limits on Y p and D/H strongly depend on fib 
(as well as on geometry and fy), ranging between about 
3 x 10~ 10 and 10~ 8 for fl b h 2 = 0.025 and 3 x 10~ 9 to 
3 x 10~ 8 for flbh 2 — 0.051, This is typically well in excess 
of the claimed primordial 7 Li/H w 1.7 x 10~ 10 as de- 
rived from observations of lithium abundances in metal- 
poor halo stars belonging to the Spite plateau. Though 
7 Li may in principle be depleted in these stars, there 
are strong arguments against this possibility, such as the 
claimed absence of intrinsic dispersion of 7 Li abundances 
in stars belonging to the Spite plateau. Further, the more 
fragile 6 Li isotope which should have been astrated as 
well, is by now observed in a few of these stars |34j . Re- 
cently, a primordial 7 Li/H even as low as 1.2 x lCp^ has 
been claimed [ [35| , which results from correcting 7 Li abun- 
dances for galactic cosmic ray production of this isotope. 
Nevertheless, even an SBBN scenario (at flbh 2 = 0.02) 
yields 7 Li/H w 3.8 x 10~ 10 , in excess of the primordial 
7 Li/H determination, when one assumes D/H 3 x 10~ D 
as favored by the QAS data. One thus would have to 
resort to a small amount of 7 Li depletion and/or system- 
atic errors due to, for example, the use of inappropriate 
stellar atmospheric models, to reconcile these values. If 
one demands the IBBN yield of 7 Li to be at, or below, 
the quoted SBBN reference value, one finds that the flbh 2 
should be below 0.025, precluding a substantial increase 
of SI;, over that inferred from SBBN. Only if one were 
to relax the 7 Li/H limit to about (seemingly unreason- 
able) 10 -9 could flbh 2 in IBBN scenarios be consistent 
with the currently preferred flbh 2 from CMBR anisotropy 
measurements ~ 0.03. 

We have tested if these conclusions could be changed 
due to existing reaction rate uncertainties quoted by 
the NACRE collaboration [0. We have changed the 
following reactions to their quoted limits: D(p,7) 3 He 
(lower limit), 3 H(a,7) 7 Li (lower limit), 3 He(a,7) 7 Be 
(lower limit), 7 Li(p,a) 4 He (upper limit), and D(a,7) 6 Li 
(upper limit). These changes have been "designed" to 
minimize Li production and maximize 6 Li production. 
With these modified rates we have performed two calcu- 



lations (a) spherical shell, flbh 2 = 0.038, /■ 



1/3 



0.025, 



iioo = 724 m, and (b) spherical condensed, flbh 2 = 0.038, 
32 m, where the length scales have 



f v /3 = 0.0125, ; 100 



been chosen close to the "optimum" distance for mini- 
mum 7 Li production. This has lead to a 7 Li/H yield of 



3 




io- 
io- 

10 
10" 

io- 
io- 
io- 

10" 

io- 

10" 
10" 



= 1 1 1 II II 1 1 1 1 1 1 1 1 

r 


i 1111111= 


— i i i 


6 Li/Hl 
i i ■*= 




10° 



10 1 10 2 
J ioo [ m ] 



10 3 



FIG. 1. Abundance yields of D, 4 He, 6 Li, 7 Li, 9 Be, and n B in IBBN scenarios as a function of the inhomogeneity length scale 
hoo (given in meters at T — 100 MeV). The calculation assumes spherical condensed inhomogeneities with high-density volume 
filling fraction fv = 0.125 3 and density contrast between high- and low- density regions of R — 200/ fv (see text for details). 
Except for the 4 He abundance which is given as mass fraction, Y p , all abundances are given as number fractions relative to 
hydrogen as indicated in the panels. The solid, dotted, short-dashed, and long-dashed lines refer to results for Qbh 2 — 0.012, 
0.025, 0.038, and 0.051, respectively. 
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FIG. 4. As Figure 1, but for spherical shell geometry and fv = 0.8 
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1.45 x 10 -9 (compared to 2 x 10 -9 when the central values 
of the reaction rates are used) in case (a) and 1.02 x 1CP 8 
(compared to 1.4 x 10 -8 ) in case (b), illustrating that the 
uncertainty in the predicted 7 Li remains within bounds. 
Nevertheless, a large uncertainty exists in the prediction 
for the 6 Li abundance: we have found a factor ~ 3 and 
4 increase in case (a) and (b), respectively. 

It should be interesting to explore if the abundance 
yields of 6 Li, 9 Be, and U B in IBBN scenarios may be 
as large as the abundances of these isotopes observed 
in the most metal-poor stars to date. Such a com- 
parison could yield, in principle, independent confirma- 
tion/rejection of IBBN scenarios. In the parameter space 
where IBBN yields are consistent with observationally 
inferred limits on He and D/H, we find production of 
6 Li/H ~ 10~ 14 — 7 x 10~ 13 (including the large reac- 
tion rate uncertainty in the D(a,7) 6 Li rate), implying 
a typical "maximum" enhancement factor for this iso- 
tope of about 10-30 compared to a SBBN scenario at 
the same fi^. For the 9 Be and 11 B isotopes one finds 
ranges of 9 Be/H - 10~ 18 - afew x 10~ 15 , and n B/H 

- 10" 16 - 10~ 13 . Typical IBBN yields of these iso- 
topes seem therefore still much below the observed 6 Li/H 

- 7 x 10" 12 @, 9 Be/H - 5 x 10" 14 HQ, and n B/H 
~ 10~ 12 |57j in the lowest metallicity stars to date where 
such observations have been performed. These observa- 
tions are thus inconclusive with regards to a validation of 
IBBN scenarios. Though it is not easy to completely rule 
out the possibility that there indeed exist very specific 
initial conditions for the baryon inhomogeneities which 
yield primordial production of 6 Li, 9 Be, and n B in abun- 
dance as high as currently observed in the lowest metal- 
licity stars, it seems clear that this is not the typical case. 

In summary, we have performed numerical simulations 
of BBN in the presence of an inhomogeneous baryon dis- 
tribution for wide ranges of the parameters describing 
the inhomogeneities and for a few representive baryon- 
to-photon ratios. Our choice of initial conditions is lim- 
ited to scenarios where essentially all baryons are within 
overdense pockets and the remainder of the volume is ini- 
tially void of baryons. We found that such scenarios may 
be consistent with observational limits on the primordial 
4 He and D abundance for fif,ft, 2 as large as ~ 0.05, how- 
ever, they result in significantly overabundant production 
of 7 Li with respect to the 7 Li/H ratio as observed in stars 
belonging to the Spite plateau. We note here that similar 
conclusions have been recently drawn by Ref . |38) . Typi- 
cal production of 6 Li, 9 Be, and n B in such scenarios are 
found to be still below the abundances of these isotopes 
observed in the most metal-poor stars to date. Unless 
'Li in stars on the Spite plateau has been significantly 
astrated, which seems unlikely, IBBN scenarios thus do 
not allow for a significant increase of fif, over that inferred 
from a SBBN scenario. 

We wish to acknowledge several useful discussions with 
In-Saeng Suh and Naoki Yoshida. 



APPENDIX A: REEVALUATION OF BARYON 
DIFFUSION CONSTANTS 

In this appendix we summarize the baryon diffusion 
constants which we used in our inhomogeneous Big Bang 
nucleosynthesis calculations. Some of these diffusion con- 
stants have been reevaluated. Such a reevaluation seemed 
necessary not only since prior work on the subject [p]j4C|l 
yielded partially conflicting results (e.g., the proton dif- 
fusion constant due to proton-electron scattering as com- 
puted by Applegate, Hogan, & Scherrer and Banerjee & 
Chitre), but also due to improvement on approximations, 
such as an energy independent neutron-proton cross sec- 
tion. Furthermore, we correct for the electron diffu- 
sion constant due to electron-photon scattering as given 
in |19|. Rather than going over the partially lengthy de- 
tails of the calculations we performed, we will state our 
results, outline by what procedure we obtained them, and 
highlight the differences to prior evaluations. 

Banerjee & Chitre j4(j] (hereafter; BC) computed 
diffusion constants by using the first-order Chapman- 
Enskog approximation for arbitrarily relativistic particles 
as thoroughly discussed in the monograph by de Groot, 
van Leeuwen, & van Weert jy] (hereafter; GLW). The 
master equation given in BC for the computation of diffu- 
sion constants (i.e., Eqs. (1), (3), and (4) in BC) are not 
directly evident from GLW but involve a fairly detailed 
computation. We have therefore redone the calculation 
of this master equation and arrive at the same result |39| ] 
as BC. Note that the first-order Chapman-Enskog ap- 
proximation is typically accurate to within 20-30 % |ll| . 

(a) neutron-electron scattering 

At higher temperatures (T > 50— 100 keV), and when 
the local baryon-to-photon ratio (77) is not too large, 
the diffusion of neutrons is limited by magnetic moment 
scattering off electrons and positrons. Using the mas- 
ter equations of BC, under the assumption of an energy- 
independent cross section, and to lowest non-trivial order 
in the small quantities m e /mN and T/m^v, where m e , 
to at, and T are electron mass, nucleon mass, and tem- 
perature, respectively, but for arbitrary T/m e , we find 



D r . 



3 /tt 1 1 K 2 {z e ) 
8 V 2« e± < e i/2^ 5/2 ( Ze ) 



(Al) 



in agreement with BC. In this expression z e — m e /T, the 
quantity n e ± is the total number density of electrons and 
positrons, the transport cross section is 

oL = Jdn^(l-cos9), (A2) 

to be evaluated in the center-of-mass system, and K n are 
modified Bessel functions of the second kind and of n'th 
order, i.e. 

1 poo 

K 2 (z) — n I k 2 exp (-V^ + z 2 ) dk , (A3) 
z 2 Jo 
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and 



K 5/2 (z) 



7T _ 7 1 3 3 

— e z 1 + - + 



(A4) 



The expression Eq. (Al) does agree with that derived 
by Applegate, Hogan, & Scherrer B (hereafter; AHS) 
via considering the drag force exerted by on neu- 
trons and using the Einstein relation. As noted by 
Ref. 0), both derivations of the diffusion constant ap- 
proximate the fermionic occupation number by a rel- 
ativistic Maxwellian, i.e. / = (exp(E/T) + w 
exp(—E/T), which nevertheless should only result in a 
small error. Using a l ne = 37r (cm/rajv) 2 , with a the fine 
structure constant and K — —1.91, we may give a numer- 
ical value for the diffusion constant 



D„ 



1.87 x 10 



4 m V T \ 



1/2 



K 2 (z e ) 



\MeVj K±/MeV 3 ) K 5/2 (z e ) 



(A5) 



where n e ± is given in natural units, i.e. h = c = 1. The 
correct density n e ± is easily obtained from the BBN code. 

Recently Suh & Mathews (i^Q] have considered finite- 
temperature effects on the neutron diffusion constant. 
They find a transport cross section due to neutron- 
electron (positron) magnetic moment scattering which 
significantly increases over a l ne — 3Tr(aK/m^) 2 at low 
T < 0.5 MeV. Since finite-temperature effects should 
vanish in the limit T — > their result is fairly surpris- 
ing. We have therefore reevaluated the neutron-electron 
cross section by using the result given in Ref. [Q and 
confirm that neutron-electron scattering is independent 
of energy for electron energies much below the nucleon 
mass. Similarly, the authors claim a significant increase 
of D ne at high T j| . Within the context of their analysis, 
such an increase could only occur due to a change in the 
electron mass and/or the transport cross section. Never- 
theless, according to their own analysis both quantities 
don't seem to deviate much from their zero-temperature 
limits for temperatures below T < 3 — 5 MeV. In light of 
these inconsistencies we therefore prefer to use the stan- 
dard AHS and BC results in our calculations. 

(b) neutron-proton scattering 

Neutron-proton nuclear scattering limits diffusion of 
neutrons at lower temperatures and/or high rj. At the low 
energies relevant for BBN the scattering cross section is 
dominated by scattering-angle independent s-wave scat- 
tering (zero angular momentum) resulting in a transport 
cross section which equals the total cross section, 



3ira 2 t 



' np 



(a s kf + (1 - \r s a s k 2 ) 2 (a t k) 2 + (1 - \r t a t k 2 ) 2 

(A6) 

In this expression k is the nucleon wave vector in the 
center-of-mass (!) system [Es), an d the parameters 



a s = -23.71 fm, a t = 5.432 fm, r s = 2.73 fm, and 
r t = 1.749 fm. Only at very low temperatures the above 
cross section is approximately independent of energy. In 
that case, one may derive (ij 



3V^ 1 / T 

8 n p a np \m N 



1/2 



(A7) 



for the diffusion constant, which is a factor of two smaller 
than the result given in BC. Nevertheless, the energy 
dependence of the cross section becomes large at T > 
50keV and should be taken into account. This may be 
done properly by evaluating the diffusion constant via 
the master equ ations in BC with the appropriate cross 
section Eq.(A6). Following this procedure we derived the 
cross section 



,- 5 mV T \ 



1/2 



1 



D nv = 2.82 x If) -0 — I - — ; 

p s \MeV J (n p /MeV 3 ) 

1 



J(oi,6i) + 0.16/(02, 6a) 



where 



1 f°° 
I(a,b) = - dx 



oaf +(1- bx/2) 2 ' 



(A8) 



(A9) 



are integrals to be evaluated numerically. Here the pa- 
rameters a and b are given by a\ = 13.59 (T/MeV), 
6i = -1.56 (T/MeV), a 2 = 0.71 (T/MeV), and b 2 = 
0.23(T/MeV). We note that the integral I(a, b) con- 
verges only slowly against its limiting value 1(0, 0) = 1 
as the temperature is decreased. In that limit, Eq. (|A8|) 



converges against Eq. (A7) with a 



' np 



37T0 2 . 



(c) proton-electron scattering 

The diffusion of protons in IBBN scenarios is only sig- 
nificant at lower temperatures. As long as Debeye screen- 
ing of proton charge in the plasma is effective, protons 
may diffuse independently of the additional (net) elec- 
trons required by charge neutrality (cf. electron-photon 
scattering). In that case, proton diffusion is limited by 
Coulomb scattering on e ± . Both, AHS and BC compute 
the proton diffusion constant due to Coulomb scattering 
D pe . However, their results differ by as much as a factor 
of eight at low temperatures. We have therefore recom- 
puted Dpe by using the master equations given in BC and 
employing the Mott scattering cross section. To lowest 
order in m e /mN, and accurate to first order in T/m e , we 
obtain 



T 2 



4V27T a 2 n e ± 



T 



1/2 i I 15 T 
1 ~r 8 m c 

A + 2-£.(A- 1) ' 



(A10) 



where A is the well-known Coulomb logarithm with Asj 
ln(T 2 m e /2nan e ±) 1 / 2 . In the limit T/m e -> Eq. flAlCj) 



reproduces the result of AHS. We conclude that D pe as 
calculated by BC is a factor of eight to small at low T. 
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Deviations between Eq. (A1C) and the result of AHS to 
first order in T/m e are due to AHS approximating the 
electron energy in the Mott scattering cross section by 
m e . Eq. ( A10 ) yields for the numerical value of the pro- 
ton diffusion constant 



D pe = 9.29 x 10" 



1 



rp \ 5/2 

MeV/ (n e ± /MeV 3 ) 



(A/5) + ± ((A/5) -1/5) 



(AH) 



(d) electron-photon scattering 

With the decrease of temperature thermally produced 
electron-positron pairs become rare and Debeye screen- 
ing of nuclear charge becomes inefficient. In this limit, 
electric forces which would rapidly be build up if the 
proton- and (net) electron- distributions differed, prevent 
the independent diffusion of these two species |L9| . One 
may show, by evaluating the electric fields which would 
be present due to differing proton- and electron- distribu- 
tions in the presence of Debeye screening, and by compar- 
ison of the resulting proton flux due to the electric fields 
with the flux of protons due to diffusion, that protons 
may only diffuse independently, when n e ± 3> n e - — n e + . 
When this is not the case, electrons and protons diffuse 
together by ambipolar diffusion, with the effective diffu- 
sion constant given by twice that of the larger of electron- 
and proton- diffusion constants jl6| . Electron diffusion is 
rendered fairly inefficient due to Thomson scattering of 
electrons on the cosmic background photons. The diffu- 
sion constant may be computed by considering the drag 
force on an electron due to a photon blackbody |f7j 



/dra 



3l5 



a T hT v , 



(A12) 



with <7xh ~ 6.65 x 10~ 29 m 2 the Thomson cross section 
and v the velocity of the electron in the cosmic back- 
ground photon rest frame. Note that Eq. (A12) is given 
in natural units. Using Eq. ( A12 ) together with the Ein- 
stein relation D = Tb, where the mobility b is defined 
as the proportionality constant between the terminal ve- 
locity v which a particle reaches in a plasma when an 
external force / is applied, i.e. v = bf , one finds for the 
effective proton diffusion constant 



2D P 



( T 



7.86 x 10~ z — . 

s VMeV 



-3 



(A13) 



applicable when n e + < n e ~ . We note here that Eq. ( [A 13] ) 
is different from the simple estimate given in Ref. |p.9|. 
We stress that the neglect of D° due to Thomson scat- 
tering, typically important at low T < 40 keV, may lead 
to errors in the calculated 7 Li abundances by more than 
an order of magnitude. 
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